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Abstract 

We investigate the possibility of detecting Majorana neutrinos at the e^e~ Next Linear Colhder 
(NLC). We study the /J/^ + jets {Ij = e,fi,T) final states which are, due to leptonic number 
violation, a clear signature for intermediate Majorana neutrino contributions. Such signals (fi- 
nal leptons of the same-sign) are not possible if the heavy neutrinos have Dirac nature. The 
interactions between Majorana neutrinos and the Standard Model (SM) particles are obtained 
from an effective Lagrangian approach. As for the background, we considered the SM reaction 
e'^e~ — )• W~^W~^W~W~ , with two W's decaying into jets and two W's decaying into + z^(^'), 
producing extra light neutrinos which avoid the detection. We present our results for the total 
cross-section as a function of the neutrino mass and the center of mass energies. We also show the 
discovery region as a function of the Majorana neutrino mass and the effective coupling. 
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I. INTRODUCTION 



The standard model of particle physics (SM) only contains left-handed neutrinos, which 
makes it not possible to generate mass for them. One very important discovery in the field 
is the neutrino oscillations, which requires the neutrinos to posses a small mass (m^ > 0.01 
eV). Thus we need to go beyond the SM to solve this issue. One way to do it, is by the 
seesaw mechanism, which requires one o more right-handed neutrinos, generically z/^, with 
a mass term 

>C'"'^'' = Mpr-L^Yur + h.c. , (1) 

where L denotes the left-handed lepton doublet, Y the Yukawa coupling matrix , the 
doublet higgs boson and M the Majorana mass matrix. There are many extensions of the 
SM (Left-right symmetric model, SO(IO), E6, ...) with extra right-handed neutrinos which 
are singlets of the SM gauge group and for which the Majorana mass terms are naturally 
allowed [ij. 

Upon diagonalization of the mass term we obtain, besides the light neutrinos, heavy 
Majorana neutrinos (N), which allow for lepton number violation. 
By solving the eigenvalue problem, we obtain the masses 

TTLj^ = mxtM^^m^, with m^) = Y—= , (2) 

V2 

and the mixing angle Uin ~ mxi/M. In typical seesaw scenarios, the Dirac mass term m-o 
are expected to be around the electroweak scale (then Y ~ 0(1) in Eq.(l2])), whilst the 
Majorana mass M being a singlet under the SM gauge group may be very large, close to 
the Grand Unification Scale. Thus, the seesaw mechanism can explain the smallness of the 
observed light neutrino masses {uij^ ~ 0.01 eV) and clearly leads to the decoupling of A^. 
Even a different choice in which M ~ 100 GeV and mv ~ 0.1 me, keeping rriy ~ 0.01 eV, 
implies a vanishing mixing Uin ~ 10^^ ^. 

This mixing weighs the coupling of with the standard model particles and in particular 
with the charged leptons through the V — A interaction: 



Cw = --^UinN'YPlIW; + h.c (3) 

This effect is so weak that the observation of Lepton Number Violation (LNV) must indicate 
new physics beyond the minimal seesaw mechanism, as was indicated in Ref.|2j. In view of 



the above discussion we consider, in a model independent way, the effective interactions of 
the Majorana neutrino A^, of mass lower than 1 TeV and negligible mixing to vl . 

In the case that heavy neutrinos (A^) exist, the present and future experiments will be 
capable of determining their nature. The production of Majorana neutrinos via e"'"e~, e~7, 



12|. 



77 and hadronic collision have been extensively investigated in the past 

In this work we study the possibility for the e^e~ next linear collider (NLC) to produce 
clear signatures of Majorana neutrinos in the context of interactions coming from an effective 
lagrangian approach. We study the reaction e+e^ — )■ /J/^ + jets {Ij = e,/i, r), which is 
divided into two subprocesses depicted in FigH] and Fig|2J In the first case we produce two 
Majorana neutrinos (A^) which will decay into one charged lepton and jets (A^ — )■ / + jets). 
In the second case, which is a three body reaction, we consider single neutrino production 
decaying in the same way as before, and a W decaying into two jets {W — )■ jets). We 
have not considered the pure lepton channels since they involve light neutrinos which escape 
detection, in which case the Majorana nature of the heavy neutrinos would have no effect 
on the signal (we should be able to know whether the final state contains neutrinos or 
antineutrinos) . For the decay of the Majorana neutrinos we have calculated the branching 
ratios of the most important channels. It is possible to parameterize the effects of new 





FIG. 1: Diagram contributing to the production of two Majorana neutrinos. 



physics beyond the standard model by a series of effective operators O constructed with 
the standard model and the Majorana neutrino fields and respecting the Standard Model 
SU{2)l (8) U{1)y gauge symmetry These effective operators represent the low-energy 
limit of an unknown theory. Their effects are suppressed by inverse powers of the new physics 
scale A for which we take the value A = 1 TeV. Here we consider the effect of dimension 6 
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operators which are the dominant. 

The total lagrangian is organized as follows: 



£ = £ 



SM 



(4) 



where J7 and i labels the operators and families respectively. For the considered operators 
we follow Ref ^| starting with a rather general effective lagrangian density for the interaction 
of a Majorana neutrino with leptons and quarks. All the operators which we list here are 
of the dimension 6 and could be generated at tree level in the unknown fundamental high 
energy theory: 



and for the baryon-number conserving 4-fermion contact terms we have: 



e,; 



(5) 



(6) 
(7) 
(8) 
(9) 



where e,, Wj, di and Lj, Qi denote the right handed SU{2) singlet and the left-handed SU{2) 
doublets, respectively. 

The operators listed above contribute to the effective lagrangian 

A^ e" A^ 





FIG. 2: Diagram contributing to the production of a single Majorana neutrino. 
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asli'^L,iNReL,ieR,i - eL,iNRi'L,ieR,i) + a^sl{'^L,iUR,iNvL,i + dL,iUR,iNeL,i) + 
Oisl{vL,iNRdL,idR^i - eL,iNRUL,idR^i) + agl{uL,iNReL,idR^i - dL,iNRiJL,idR,i) + 
asli^Ri^L,i^L,iNR + iVKeL,iei,,iVR) + /i.e.} (10) 

where the sum over i is understood and the constants a^j are associated to specific operators 
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(11) 

We calculate the cross-section for the production of the Majorana neutrino according to the 
process shown in FiglH valid for the kinematic range m^r < y/s/2 



cr 



NN in- , , 



647rA4 



where 



+ (12) 



, ^ (13) 

and Cl = —1/2 + xvi/, = xw, with x^y = sin^6'pi/. We also study the single Majorana 
neutrino production. The corresponding diagrams are shown in Figj2] and the result for the 
square amplitud is 

I-MP = |8a^^.[2(fc ■ p){2{k .q){l.q)-{k- l)q') + M^(2(/ ■q){p-q)-{l- p)q')] 

+ 4(4?)'^[2(^ ■ ■ q){p, ■ q) - {k ■ p,)q') 

w 

+ Ml{2{p-q){p^-q)-{p-p^)q^)]] (14) 

where q = p + k and pi, p2, I, p and k are the 4-momenta of the electron, the positron, the 
Majorana neutrino A^, the charged lepton and the W boson respectively. The W propagator 
is Uw = ml,/[{pi - If - m^]. 
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FIG. 3: In the left panel the Branching ratios for the Majorana neutrino decay with = 0. In 

the right panel there is a comparison with the same-coupling case for the N ^ + jets decay. 
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FIG. 4: The cross-section for the process e~^e~ — )• l^l^ + jets. The left panel correspond to the 
zero contribution of the operators related with Ovijp {aQvpp =0). In the right panel the curves 
labeled (a), (b), (c) and (d) correspond to energies of center of mass y/s = CSTeF, ^/s = ^.KTeV 

(i) 

and different values of the constant aj . In both panels the horizontal solid line represent the 
values of the SM background. 



The cross section is obtained by integrating the phase space in the usual way, using the 
numerical routine RAMBO jl^ . 

The total cross section is the combination of the two processes mentioned above in the 
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approximated expression: 

^ie+e-^i±i±jets) ^ ^ ( ^ a(^^^"^^^)5r(iV ^ /+ + jets)Br{N ^ /+ + jets)e(v^/2 - m^v) 

+ ^ ltjets)Br{W jets)e(v^ - m^v) j (15) 

The factor two in front of Eq.( |T5ll takes into account the possible charges of the final leptons. 

As we shall see later, some of the considered operators contribute to the neutrinoless 
double beta decay (Oz/^/3-decay) and may be strongly constrained. In these conditions we 
analyze the case where the operators contributing to the Ou^jj-decaj vanishes whist the rest 
are non-zero and contribute with similar strength. For completion we also analyze the case 
where all of the operators contribute with similar strength. The branching ratios shown in 
the left panel of Fig. ([3]) (the expressions are collected in the Appendix) correspond to the 
case with non contribution of the Oi/^^-decay related operators {aQyjjp = 0). In the right 
panel of the same figure we show for comparison the branching ratio for N + jets with 
(^ovisp = and aQ^pp 7^ 0. As we can see there are not significant differences. 

In FigJUwe show the results for the cross section combining the processes shown in FigsH] 
and [2] with the VT-boson decaying into hadrons and the Majorana neutrino decaying 
according to the Branching Ratios shown in the Appendix. We show the result as a function 
of the Majorana neutrino mass m^r and center of mass energies of y/s = 0.5 TeV and 
^/s = 0.8 TeV. We have considered ^/s < A in order to ensure the validity of the effective 
lagrangian approach. The cross section is calculated for different values of the constants 
. In the left panel of the FigJH we have shown the cross section for the case in which the 



a 



constants related with the operators contributing to Oz/^j^-decay are considered zero. These 
operators are O]^^^, 0\^^^, ^quNL^ ^lnqcI ^qnlcLi "^i^^ ^e discussed in section lllli 
In the right panel we plot the cross section where the operators that contribute are the 
4-fermion operators (solid line) or the operators involving bosons (dashed line). In both 
panels, the non-zero coupling constants take the value one. As we can see, the 4-fermion 
contribution is the most important. In both panels we show with a horizontal solid line the 
value of the SM background as we will be explained later in the text. 

The final leptons can be either of e^, /i^ or since this is allowed by the interaction 
lagrangian (Eq JTOl) . All of these possible final states are clear signals for intermediary Ma- 
jorana neutrinos, thus we sum the cross section over the flavors of the final leptons. The 



partial width of was determined at tree level considering the dominant decay modes 
— 7- Z + 2jets, N-^l + tb, N-^l + leptons, N ^ u + H and N ^ u + quarks coming 
from the higgs, the charged VT-boson and the 4-fermion effective interactions. We present 
in the Appendix the differential partial-width for its dominant decay channels, where the 
contributing effective operators are identified by the indicated labels in the couplings. 



II. THE STANDARD MODEL BACKGROUND AND THE DISCOVERY RE- 
GION 



The considered signal is strictly forbidden in the Standard Model. The SM background, 
which can be confused with the studied reaction, will always involve additional light neu- 
trinos. The dominant SM process arises from the resonant production of four bosons: 
e+e" W'^W^W'W', the decay of two W's into leptons + and the 

other two into jets, W jets. We calculated the cross section for these processes using 



the package COMPHEP 19l-|2l| and we multiplied it by the corresponding branching ratios 



{BR[W /z/])2 ~ 0.011 and the {Bt[W -> 2jets])^ ~ 0.46 and by the factor 18 to take 
into account the different combinations of the same sign charged final leptons: l^l^ with 
/ = e, /i, r. The calculated values are 5.0 10~^ pb and 1.1 10^^ pb for ^/s = 0.8 TeV and 0.5 
TeV respectively. In Table H] we compare the values of the signal, for different values of the 
Majorana neutrino mass, with the SM background. In FigJUwe show along, with the signal 
cross-section, the corresponding background levels as horizontal lines for ^/s = 0.5TeV and 
= 0.8TeV. 

In order to investigate the capability of the studied process to discover effects of Majorana 
neutrinos, we study the region (discovery region) where the signal can be separated from 
the background with a statistical significance of 5a. It is done by defining the quantity S 

where L is the luminosity and the numerator represents the discrepancy between the signal 
and the SM background {(Tb)- In FigElwe show the discovery region (above the solid curves) 
where S > 5 {5 a statistical significance) for a luminosity L = 100 fb~^. 

For completion we have also considered, although in an approximated way, the bounds 
on the operators which come from the Oi/^/j-decay and from LEP and low energy data. The 



8 




N 



(1) 



(2) 



FIG. 5: Contribution to T-L in Eq. (j22p . In the diagram (1) the sohd dot represent the operator 
'^Ne<p and in the diagram (2) the 4-fermion operators O]^^^, OquTVL- ^^LNQd and O^^j^^ 



former will be considered in the next section and the latter in the following. 

The heavy Majorana neutrino couples to the three flavors families with couplings de- 
pendent on the scale A and the constants a^j . It is possible to relate this coupling with 
the mixing between light and heavy neutrinos {11^^, U^m, Utn) for which the experimental 
bounds, obtained from LEP and low energy data, have been put in 15l-ll8|. This relation 
was found in comparing the operator O]ve0 with the strength of the V-A interaction 
(Eq.dn])). It is Ui^M ~ (a[^/2)(t>^/A^) where v corresponds to the vacuum expectation value: 
V = 250GeV. In order to keep the analysis as simple as possible we consider that the same 
bound applies for all the couplings a^j , generically a. 



In our case, for one heavy Majorana neutrino, and following 



15|, we have: 



^iv — UinUi 



UN 



(17) 



where the allowed values for the parameter are 22| : 



fiee < 0.0054, Vt^,^, < 0.0096, < 0.016 



(18) 



For the Lepton-Flavour- Violating process (LFV), e.g. /i — )■ ej, fi — )• eee and r — )■ eee, which 
are induced by the quantum effect of the heavy neutrinos, we have |23| : 



^en\ < 0.0001, \ner\ < 0.02, < 0.02 



(19) 



These bounds can be translated to the constants a considering, in a simplified way, that all 
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TABLE I: Comparison between the signal and background, 
the operators satisfy the same constraint 

For A = 1 TeV we have: 

a < 0.32 

This value is shown in both panels of Figj6] with a horizontal dot-dashed line. 



2 \ 2 
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(20) 
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FIG. 6: Discovery region above the sohd curve at 5 o" (this work) and below the dot-dashed line 
for LVF bound (Left panel) or the dotted curve for Oz^^^-decay bound (Right panel). 
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III. NEUTRINOLESS DOUBLE BETA DECAY BOUNDS 



In order to take into account the bounds imposed by the Oz/^^-decay experiment on some 
of the couphng constants , we consider, in a general way, the following effective interaction 
Hamiltonian: 

H = Geff uTd eTN + h.c. (22) 



where 
in 24, 



represents a general Lorentz-Dirac structure. Following the developments presented 



25| and using the most stringent limits on the lifetime for neutrinoless double beta 



decay To,y^^ > 1.9 x 10^^ yr obtained by the Heidelberg- Moscow Collaboration [26|, we have 
obtained the following bounds for Ge// 

The lowest order contribution to Oz/^^-decay from the considered effective operators comes 
from those that involve the W field and the 4-fermion operators with quarks u, d, the lepton 
e and the Majorana neutrino A^: 

^Nerf) ' ^duNe ' ^QuNL ' ^LNQd ' ^QNLd (24) 

The contribution of these operators to the effective Hamiltonian in eq. fl22|) is shown in FigJS] 

For the coupling constant associated with each operator we use the generic name 
that is to say 

_ (1) _ (1) _ (1) _ (1) _ (1) 

In order to estimate the bounds on the different a^j we consider two different situations. 
First, we suppose that the contribution of all the operators involved in the Oi^^/j-decay 
adds constructively. In this case we expect strong limits on the couplings, then we may 
assume them to be negligible. Thus, we consider that all of the constants associated with 
the operators in Eq. fl2^ vanish (aoi//3/3 = 0) and that the other constants, which are not 
bounded by neutrinoless double beta decay, are non-zero and have similar magnitude. This 
situation is shown in the left panel of the Fig.(l6l). 

Second, we consider the individual contributions of each operator as acting alone. In this 
case it is obvious to relate the coupling with the Ge// in Eq. (l22|) 



Ge// - (2bj 
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Thus we can translate the hmit which came from Ge// to aoup/B- For, A = ITeV, it is 

Taking a conservative point of view, in the right panel of Figj6l we present this bound 
considering that it is the same for all the constants (generically a) and show it with the 
dotted curve. The solid curve, which is the contribution of this work, represent the lower 
limit for the discovery region. In the right panel, it was calculated considering that all the 
constants have similar magnitude. On the other hand, in the left panel it was calculated 
considering aoi//3^ = 0. We also show the bound from Lepton Flavors Violating process with 
the dot-dashed line in the same figure. 

Summarizing, we calculated the cross-section for the process e'^e'^ — IJl^ + jets where 
/i, I2 and I3 are light leptons {e,fi,T) respectively. We show the total unpolarized cross- 
section using the calculated Branching ratios for different values oi tun and the coupling 
a^j . We showed the discovery regions at ha statistical significance combining with the Qvpp 
and the LFV bounds. We found that it will be possible to discover Majorana neutrinos with 
masses lower than 250 GeV and 400 GeV at e+e^ colliders with center of mass energy of 
0.5 TeV and 0.8 TeV respectively. 
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IV. APPENDIX 

We present here the partial decay widths of a heavy Majorana neutrino N for its dominant 
decay channels. They were calculated using the effective interactions shown above in the 
text. 
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with < a; < 1 - % 



{N leptons) 

dx 



L v=i,3 




with < a; < 1 



dx 



2567r3 V A 7 




277^(1 - z^) 



with < X < 1 



where x = ^p^^^^^Jtun, and y = {rrit/mNf, = {m^/mNf, z = {rriN/mwf, w = 
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